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We report the first observation of e^e^ T(15')7r''"7r~ , T(25')7r^7r~, and first evidence for 
e+e" r{3S)n+n~, r{lS)K+K', near the peak of the T{5S) resonance at ^ ~ 10.87 GeV. 
The results are based on a data sample of 21.7 fb~^ collected with the Belle detector at the KEKB 
e"'"e~ collider. Attributing the signals to the T(5S) resonance, the partial widths r(T(5S') 
T(lS)7r+7r-) = 0.59±0.04(stat)±0.09(syst) MeV and r(T(5S') T{2S)n+TT-) = 0.85±0.07(stat)± 
0.16(syst) MeV are obtained from the observed cross sections. These values exceed by more than 
two orders of magnitude the previously measured partial widths for dipion transitions between lower 
T resonances. 

PACS numbers: 13.25.Gv, 14.40.Gx 



Heavy quarkonia provide a unique nonrelativistic sys- 
tem in which low energy QCD may be illuminated 
through their energy levels, widths, and transition am- 
plitudes. Dipion transitions between ^ and T levels be- 
low the open flavor thresholds have been successfully 
described in terms of multipole moments of the QCD 
field [l|. The first measurements above the open beauty 
threshold, namely of T(45) ^ T(15)7r+7r- are 
consistent with this picture fH"]. (The T(4S') is the third 
radial excitation of the J^*-^ = 1 state T(IS').) 

The spectroscopy above open flavor threshold is com- 
plex, however, as there is no positronium analogue. The 
recent discovery of a broad 1^^ state, the y(4260), 
decaying with an unexpectedly large partial width to 
J/^Tr'^TT~ 6], has brought new challenges to the inter- 
pretation of its composition, with "hybrid" ccg (where g 
is a gluon) and ccqq (where qq is a color-octet light quark 
pair) four quark state as possibilities. The observation of 
a bottomonium counterpart to F(4260), which we shall 
refer to as Yf, [7] , could shed further light on the structure 
of such particles. The expected mass is above the T(4S'). 
It has been suggested that a Yf, with lower mass can be 
searched for by radiative return from the T(55), and one 
with higher mass through an anomalous rate of T(ri5)7r7r 
events scaling from T(4S') T(lS')7r7r, one expects 
T(55') T(lS')7r7r to have branching fraction ~ 10^^. 

Here we report the first observation of T(lS')7r+7r^ 
and T(25)7r+7r^ final states, as well as evidence for 
T(35')7r+7r- and T{IS)K+K- in a 21.7 fb'^ data sam- 
ple collected near the peak of the T{5S) resonance with 
the Belle detector at the KEKB e+e~ energy-asymmetric 
collider [Ij. The rates for T(15)7r+7r- and T(25')7r+7r- 
are much larger than the expectations from scaling the 
comparable T(4S') decays to the T(55). Since only one 
center-of-mass (CM) energy is used, one does not know 



whether these enhancements are an effect of the T(55') 
itself, or due to a nearby or overlapping "Yj," state. 
Throughout this Letter, we shall therefore use the no- 
tation T(10860) instead of T{5S). 

The Belle detector is a large-solid-angle magnetic spec- 
trometer, which consists of a silicon vertex detector 
(SVD), a central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), a barrel- like ar- 
rangement of time-of-fiight scintillation counters (TOF), 
and an electromagnetic calorimeter comprised of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
that provides a 1.5 T magnetic field. An iron flux-return 
located outside the coil is instrumented to detect 
mesons and to identify muons (KLM). The detector is 
described in detail elsewhere 0. 

The T(10860) T{nS)TT+7r- and T{1S)K+K~ final 
states are reconstructed using T{nS) — s- fJ-^fJ.'' decays. 
Events with exactly four well-constrained charged tracks 
and zero net charge are selected. Muon candidates are 
required to have hits in the KLM detector associated 
with the extrapolated trajectory of the charged track. 
Two muons with opposite charge are selected to form a 
T{nS) candidate. The two remaining tracks are treated 
as pion or kaon candidates. To suppress the background 
from /x+/i~7 — > /i+/i^e+e~ with photon conversion, pion 
candidates with positive electron identification are re- 
jected. Electron identification is based on associating 
the ECL shower energy to the track momentum, dE/dx 
from CDC, and the ACC response. Kaon candidates are 
required to have a kaon likelihood, estimated with infor- 
mation from the ACC, TOF, and dE/dx from the CDC, 
greater than 0.1. This requirement has an efficiency of 
98.2%. The cosine of the opening angle between the 7r+ 
and 7r~ {K^ and K^) momenta in the laboratory frame 
is required to be less than 0.95. The trigger efficiency is 
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found to be very close to 100% for these final states. To 
reject (radiative) Bhabha and ^-pair backgrounds, the 
data are required to satisfy either 9^3.^ < 175°, or 2 GeV 
< J2 ^ECL < 10 GeV, where ^max is the maximum open- 
ing angle between any charged tracks in the CM frame, 
and J2 ^ECL is the sum of the ECL clusters' energy. 

The signal candidates are identified using the kine- 
matic variable AM, defined as the difference between 
Min+n-TT+TT-) or M{p+n-K+K-) and Mifx+n') for 
pion or kaon modes. Sharp peaks are expected at 
AM = Mr{raS) - Mt(„s) for m > n. For T(10860) 
T(n5')7r+7r^ and T{1S)K^K~ , signal events should be 
concentrated at AM — y/s — M-Y(^nS)7 since a single CM 
energy is used. 

Figure [T] shows the two-dimensional scatter plot 
of M{iX^ii^) vs. AM for the data. Clear en- 
hancements are observed, especially for T( 10860) 
T(lS')7r+7r- and T(2S')7r+7r- decays. The dom- 
inant background processes, e+e^ fi'^ii^j{-^ 
e+e") and e+e" /i+/i^7r+7r^ accumulate at the 
kinematic boundary, M(/i"*'/i~7r~''7r~) = y^. The 
events with |M(/i+/i^7r+7r") - y/s\ < 150 MeV or 
\M{n+^i-K+K-) - Vi| < 150 MeV are selected. The 
fitting regions are defined by 1.25 GeV/c^ < AM < 



1.55 GeV/c2, 0.69 GeV/c^ < AM < 0.99 GeV/c^, and 
0.36 GeV/c2 < AM < 0.66 GeV/c^ for T(10860) ^ 
T(lS')7r+7r-, T(2S')7r+7r-, and T(3S')7r+7r-, respectively. 
The fitting region in AM for T(10860) T{IS)K+K- 
is the same as for the T(lS')7r+7r^ mode. The oblique 
fitting regions are selected so that the background shape 
is monotonic along each band. The background dis- 
tributions are verified using the off-resonance sample 
(recorded at 10.52 GeV) 0. 

The AM distributions for the //"""^"tt+tt^ candidates 
in the T(IS') and T(2S') fJ-~^fJ'~ mass bands are shown 
in Fig. [1 The peaks for T(10860) -> T(lS')7r+7r- 
and T(25)7r+7r- are located at AM ~ 1.41 GeV/c^ 
and ^ 0.84 GeV/c^, respectively. Two other peaks at 
AM - 0.56 GeV/c^ and - 0.89 GeV/c^ correspond to 
T(2S') T(15)7r+7r- and T(35) T(lS')7r+7r- tran- 
sitions, respectively. The absence of a peak around 1.12 
GeV/c^ corresponding to T(45) T(lS')7r+7r^ is consis- 
tent with the rates measured in Refs. [sj,!^. The structure 
just below T(3S') T{1S)tt^tt~ in the AM distribution 
is from the cascade decays T(10860) T(25)7r+7r^ with 

T(2S') ^ T(is')H ^J'^^^~]x. 

Signal yields are extracted by unbinned extended max- 
imum likelihood (ML) fits to the AM distributions. The 
likelihood for the fit is written as 




AM = M(^l^l7C7I)-M(n^l) (Gev/c-^ 
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FIG. 1: Scatter plot of M{fj,+n~) vs. AM for the data 
collected at y's ~ 10.87 GeV, for (a) fj,'^ fj,~ ■k'^ -jt~ and (b) 
fi'^ IJ,~ K'^ K" candidates. Horizontal shaded bands corre- 
spond to T(IS'), r{2S) and T(3S') (only T(IS') for (b)), 
and open boxes are the fitting regions for T( 10860) 
T{nS)n^Tr~ and T{1S)K'^ K~ . The lines indicate the kine- 
matic boundaries, M{fi^ h^-k^tv" , fj,'^ fj,~ K'^ K~) = y^. 



CiNs,Nb) 



N 



m 



Y[[N,-Ps{AM,)+Nb-Pb{AM^)] , 

i=l 

(1) 

where Ns (Nb) denotes the yield for signal (background), 
and Ps {Pb) is the signal (background) probability den- 
sity function (PDF). The signal is described by a sum of 
two Gaussians while the background is approximated by 
a linear function. The tail part of the signal PDF is pa- 
rameterized by a broad Gaussian, whose width and frac- 
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FIG. 2: The AM distributions for ti+ fi 
the (a) T(15') fi+fi' and (b) T(2S) fi+ fi' bands of 
Fig. [ija) . Vertical dashed lines show the expected AAd val- 
ues for the T{nS) T(l, 25')7r''"7r~ transitions. 
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tion (of area) are fixed from Monte Carlo (MC) simula- 
tion. For the T(10860) -> T(lS')7r+7r- and T(2S')7r+7r- 
modes, the remaining PDF parameters and yields of 
signal and background are floated in the fits. For the 
T(10860) ^ T(3S')7r+7r- and T{1S)K+K- transitions, 
where statistics are limited, the means and widths are es- 
tablished based on T(10860) T(lS')7r+7r~ events and 
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fixed in the fits. We observe 325_ 

and 20.2tli events in the T(10860) ^ T(lS')7r+7r-, 
T(2S')7r+7r-, T(3S')7r+7r-, and T{1S)K+K- channels, 
with significances of 20a, 14(t, 3.2a, and A.9a, respec- 
tively. The significance is calculated using the difference 
in likelihood values of the best fit and of a null signal hy- 
pothesis including the effect of systematic uncertainties. 
The Gaussian widths of the T(10860) -> T(lS')7r+7r- 
and T(2S')7r+7r- peaks are found to be 8.0 ±0.5 MeV/c^ 
and 7.6 ± 0.7 MeV/c^, respectively, and are consistent 
with the MC predictions. The distributions of AM with 
the fit results superimposed are shown in Fig. [3l 

The yields for T(10860) ^ T(lS')7r+7r-, T(2S')7r+7r- 
are found to be large; thus, the corresponding invariant 
masses of the tt+tt" system, ^(Tr+Tr"), and the cosine 
of the helicity angle, cos^Hei, can be examined in detail. 
The helicity angle, 0Hei, is the angle between the tt^ and 
T(10860) momenta in the tt+tt" rest frame. Figure [4] 
shows the T(10860) yields as functions of M(7r+7r~) and 
cos ^Hoi, which are extracted using ML fits to AM in bins 
of M{Tr~^TT~) or cos^hoI- The shaded histograms in the 
figure are the distributions from MC simulations using 
the model of Ref. [l[ , while the open histograms show a 
generic phase space model. As neither model agrees well 
with the observed distributions and the efficiencies are 
sensitive to both variables, the reconstruction cfhcien- 
cies for T(10860) ^ T(lS')7r+7r- and T(2S')7r+7r- are 
obtained using MC samples reweighted according to the 
measured M(7r+7r~) and cos^Hei spectra. Due to limited 
statistics, we estimate the reconstruction efficiencies for 
T(10860) ^ T(35>-+7r- and T{1S)K+K' modes using 
the model of Ref. [l[. Comparison of the M(7r+7r^) distri- 
bution obtained here with other T{nS) — > T(mS')7r+7r~ 
(m < n) decays could be irnportant for the theoretical 
interpretation of the results [l|, [3| • 

Assuming that signal events come only from the T(5S') 
resonance, the corresponding branching fractions and 
partial widths can be extracted using ratios to the T(55) 
cross section at - 10.87 GeV, 0.302 ± 0.015 nb 
The results, including the observed cross sections, are 
given in Table [J The values include the world aver- 
age branching fractions for T(nS) — » jX^ jJT decays, and 
the total width of the T(55') [10|. The measured partial 
widths, of order 0.6-0.8 MeV, are large compared to all 
other known transitions among T{nS) states. The par- 
tial widths for T(2S'), T(3S'), and T(4S') ^ T(lS')7r+7r- 
transitions are all at the keV level (Table |TT| . 

The systematic uncertainties for the cross sections 
are dominated by the T(nS) fJ'^ l-J'^ branching frac- 
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FIG. 3: The AM distributions for (a) T(lS')7r+7r-, (b) 
T(2S')7r+7r~, (c) T(3S')7r+7r-, and (d) T{1S)K+K- with the 
fit results superimposed. The dashed curves show the back- 
ground components in the fits. 
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FIG. 4: The T(10860) yields as functions of M(7r+7r") and 
cos^Hci for (a,c) T(15')7r^7r~ and (b,d) T(2S')7r"'"7r~ transi- 
tions. The shaded (open) histogram are from MC simula- 
tions using the model of Ref. jlj| (phase space model). The 
numerical yields are given in the appendix. 



tions, MC reconstruction efficiencies, and PDF param- 
eterization for the fits. Uncertainties of 2.0%, 8.8%, 
and 9.6% for the T(IS'), T(2S'), and T(3S') ^ 
branching fractions are included, respectively. For the 
T(lS')7r+7r~ and T(2S')7r"'"7r~ modes, the reconstruction 
efficiencies are obtained from MC simulations using the 
observed M(7r+7r^) and cos6'hc1 distributions as inputs. 
The uncertainties associated with these distributions give 
rise to 4.4% and 6.8% errors for the T(lS')7r+7r~ and 
T(2S')7r+7r^ MC efficiencies, respectively. For the other 
two channels, we try as input the models of Ref. [H 
and phase space model; the corresponding variations 
in acceptance are included as systematic uncertainties. 
A relatively large uncertainty of 13.6% arises for the 
T(10860) ^ TilS)K+K- channel, while the corre- 
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TABLE I: Signal yield (A'^s), significance (E), reconstruction efficiency, and observed cross section (a) for e'^e T(nS')7r^7r 
and T{1S)K+K- at ^ ~ 10.87 GeV. Assuming tlie T(5S') to be the sole source of the observed events, the branching fractions 
(B) and the partial widths (F) for T(55') — > T(n5')7r^7r~ and T{1S)K'^ K~ are also given. The first uncertainty is statistical, 
and the second is systematic. 



Process 




E 


Efr.(%) 


a(pb) 




r(MeV) 


T(lS')7r+7r- 




20cr 


37.4 


1.61 ±0.10 ±0.12 


0.53 ±0.03 ±0.05 


0.59 ±0.04 ±0.09 


T(2S)7r+7r- 


186 ± 15 


Ua 


18.9 


2.35 ±0.19 ±0.32 


0.78 ±0.06 ±0.11 


0.85 ±0.07 ±0.16 




10.511° 


3.2(7 


1.5 


l-44lo:45 ±0.19 


0-48lo;i5 ±0.07 


0.52l° f? ±0.10 


r{is)K+K- 


20.2l^;^ 


4.9a 


20.3 


o.i85i!5:E;4? ±0.028 


0.06llo;^i'^ ± 0.010 


0.067l^ o}5 ± 0.013 



TABLE IL The total width P 



total , 



and the partial width 



e+e- 1 J- T{lS)7r + 7i 



Most values are from Refs. [3I. [lo| . 



Process 



- total 



Y(15)7r^ 



T(25') ^ T(15')7r+7r^ 
T(3S') ^ T{lS)n+n- 
T(4S') ^ T{lS)n+n- 
T(10860) T(15')7r+7r- 110 MeV 0.31 keV 0.59 MeV 



0.032 MeV 0.612 keV 0.0060 MeV 
0.020 MeV 0.443 keV 0.0009 MeV 
20.5 MeV 0.272 keV 0.0019 MeV 



spending error for T(10860) — * T(35)7r+7r^ is small 
(3.2%) due to the limited phase space. The uncertain- 
ties from PDF parameterization are obtained either by 
replacing the signal PDF with a sum of three Gaussians, 
or by replacing the background PDF with a second-order 
polynomial. The differences between the fit results ob- 
tained with alternative PDFs and the nominal results 
are taken as the systematic uncertainty. The selection 
criteria for rejecting radiative Bhabha and /i-pair events 
are examined using the data (I4I collected at the T(3S') 
resonance. The 1.9% difference between data and MC ef- 
ficiencies for T(3S') T(lS')7r+7r~ decays is included as 
a systematic uncertainty. Other uncertainties included 
are: tracking efficiency (1% per charged track), muon 
identification (0.5% per muon candidate), electron re- 
jection for the charged pions (0.1-0.2% per pion), kaon 
identification (1.8% per kaon), trigger efficiencies (0.9- 
4.5%), and KEKB luminosity (1.4%). The uncertain- 
ties from all sources are added in quadrature. The to- 
tal systematic uncertainties are 7.5%, 13.5%, 13.1%, and 
15.3% for the T(15)7r+7r-, T(2S')7r+7r-, T(35)7r+7r-, 
and T(15)if±i4r~ channels, respectively. 

For branching fraction estimation, the error in the 
T{5S) cross section (±0.015 nb) gives a 5.0% uncertainty 
in signal normalization. For the partial widths, there is 
an additional uncertainty of 11.8% coming from using the 
total width of the T(55). 

In summary, we report the first observation of e±e^ — > 



T(15')7r+7r- and T(25)7r+ 
dence of e±e~ T(3S')7r^ 



transitions, and first evi- 
and TilS)K+K- transi- 
tions at a CM energy near the T(5S') resonance of y/s ~ 



10.87 GeV. Clear signals are observed at the expected 
CM energy, with subsequent T{nS) fJ^'^fJ-' decay. The 
measured cross sections are 1.61 ± 0.10 ± 0.12 pb, 2.35 ± 
0.19±0.32 pb, 1.44l°-^^±0.19 pb, and 0.185l|^ |^ti±0-028 
pb for e+e- T(lS')7r+7r-, T(2S')7r±7r-, T(35)7r+7r^, 
and T{1S)K+K~ transitions, respectively. The first un- 
certainty is statistical, and the second is systematic. As- 
suming the observed signal events are due solely to the 
T(55') resonance, branching fractions are measured to 
be in the range (0.48-0.78)% for T(n5')7r+7r~ channels, 
and 0.061% for the T{IS)K+K- channel. The corre- 
sponding partial widths are found to be in the range 
(0.52-0.85) MeV for T(n5')7r+7r-, and 0.067 MeV for the 
T{1S)K+K~ mode, more than two orders of magnitude 
larger than the corresponding partial widths for T(4S'), 
T(3S') or T(25') decays. The unexpectedly large partial 
widths disagree with the expectation for a pure bb state, 
unless there is a new mechanism to enhance the decay 
rate. A detailed energy scan within the T{5S) energy 
region would help to extract the resonant spectrum; a 
comparison between the yield of T(n5')7r+7r~ events and 
the total hadronic cross section may help us to under- 
stand the nature of the signal. 
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APPENDIX: YIELDS OF T(10860) ^ T(15')7r+7r- 
AND T(2S)7r+7r- TRANSITIONS AS FUNCTIONS 

OF A/(7r+7r") AND cos fed 



TABLE III: Numerical yields of T(10860) ^ T(15)7r+7r- and 
T(2S)7r^7r~ transitions as functions of M(7r"'"7r~) and cos 6uci 
which are shown in Figure [l] The uncertainties of the yields 
are statistical only. 



Bin 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 



Bin 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 



T(15')7r+7r- 
M{7v+7v-) (GeV) 



[0.25 
[0.35 
[0.45 
[0.55 
[0.65 
[0.75 
[0.85 
[0.95 
[1.05 
[1.15 
[1.25 
[1.35 



[-1.0, 
[-0.8, 
[-0.6, 
[-0.4, 
[-0.2, 
[+0.0, 
[+0.2, 
[+0.4, 
[+0.6, 
[+0.8, 



0.35) 
0.45) 
0.55) 
0.65) 
0.75) 
0.85) 
0.95) 
1.05) 
1.15) 
1.25) 
1.35) 
1.45) 



1 

25 
13 
11 
18 
28 
34 
34 
28 
60 
50 
14 



Yield 

oo+l.«4 

^•^-i.ie 

nq+5.70 
^'J_4 97 

"'-3.94 
91 +3.99 
^-■--3.28 
^^+5.12 

"-'-5.47 
8/1+6-70 
^4_5.97 
qi:+6.57 
•^^-5.86 
40+6.03 
^^-5.35 
91+8.57 

99 
00 



7.87 
+7.48 
6.81 
+4.09 
3.41 



Yield 



-0.8) 
-0.6) 
-0.4) 
-0.2) 
+0.0) 
+0.2) 
+0.4) 
+0.6) 
+0.8) 
+1.0) 



37 
37 
29 
21 
33 
23 
32 
37 
31 
31 



08-6.24 
Qg+6.90 

qr+6.07 
yO-5.41 
7^ + 5.51 
''±-4.84 
77+6.48 
' '-5.85 
pr + 5.55 
"^-4.84 
8n+6-34 
OU-5.67 
79+6.86 
'^-6.17 
Q7+6.33 
y '-5.65 
7r+6.28 
'^-5.58 



T(25')7r+7r- 
M{n+n-) (GeV) 



Yield 



[0.25, 
[0.30, 
[0.35, 
[0.40, 
[0.45, 
[0.50, 
[0.55, 
[0.60, 
[0.65, 
[0.70, 
[0.75, 
[0.80, 



0.30) 
0.35) 
0.40) 
0.45) 
0.50) 
0.55) 
0.60) 
0.65) 
0.70) 
0.75) 
0.80) 
0.85) 



cos f Hcl 



[-1.0: 
[-0.8. 

[-0.6. 
[-0.4. 
[-0.2. 
[+0.0. 
[+0.2. 
[+0.4. 
[+0.6. 
[+0.8. 



-0.8) 
-0.6) 
-0.4) 
-0.2) 
+0.0) 
+0.2) 
+0.4) 
+0.6) 
+0.8) 
+1.0) 



8 
20, 
13, 

9 
20, 
24, 
35, 
24, 
10, 

7 

4 



48^ 



0.0 

+3.51 
2.82 
qi+5.22 
■Jj^-4.52 
gg+4.38 

,37j 



3.70 
+3.73 
3.02 
CI +5.27 
^±^-4.58 
7n+5.70 
' "-4.97 
co+6.83 
^°-6.10 
7f\+5.58 
' "-4.88 
07+3.74 



3.03 
+3.35 
-2.64 
+2.58 
-1.88 



Yield 



30 
20 
18 
16 
7 
14 
6 
7 
24 
32 



■0^-5.46 
09+5.14 

■°^-4.49 
1 r+4.83 

■1^-4.14 
70+4.78 

■ '"-4.07 

.4511:?! 
90+4.53 

■^"-3.76 
r,r,+3.07 

■""-2.36 
1 c:+3.44 

40+5.65 
■^^-4.95 

9f,+6.41 
■^"-5.65 



